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The effect of different isospin-dependent cross-section on directed flow is studied for variety of 
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Abstract 



systems(for which experimental balance energies are available) using an isospin-dependent Quan- 
tum Molecular Dynamic (IQMD) model. We show that balance energies are sensitive towards 
isospin-dependent cross-sections for light systems, while nearly no effect exist for heavier nuclei. 
A reduced cross-section a = 0.9ctnn with stiff equation of state is able to explain experimental 



balance energies in most of systems. A power law behaviour is also given for the mass dependence 

o ■ 

of balance energy, which also follow N/Z dependence. 
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Q\ ■ PACS numbers: 25.70.-z, 25.75.Ld 
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I. INTRODUCTION verse in-plane flow, one has also proposed, e.g 



^ . differential [12( and elliptical flow [13j etc. 

H ' 

The heavy-ion physics branch has been In general, collective flow in heavy-ion colli- 



renewed interest very recent 
from the fusion probabilities 



y. This ranges sions is affected by both the nuclear mean 

1] to symmetry field potential and nucleon-nucleon (NN) 

energy dependence at intermediate energies cross-sections. One should also keep in the 

as well as fragmentation of colliding matter mind that reaction dynamics depends also on 

. One observable that has been used exten- the incident energy as well as on the impact 

sively for extracting information from heavy- parameter of the reaction At low in- 

ion collisions is the collective in-plane flow cident energies, reaction dynamics is domi- 



of various particles 



-111). Apart from trans- nated by the attractive nuclear mean field po- 
tential which results in deflection to negative 



Electronic address: suneel.kumar@thapar.edu 



angle. Worth mention, at these enegies the 
phenomenon fusion, and cluster radioactivity 
are dominated jl[ [3] . With increasing inci- 
dent energy, repulsive nucleon-nucleon scat- 
tering becomes important and results in re- 
duced negative flow caused by the attractive 
mean field potential. As a result, at a certain 
incident energy, called the balance energy, 
the in-plane flow vanishes result of can- 
cellation between these two competing effects 
[10-14]. The composite dependence of the 
Et, a i on the mean field and nucleon-nucleon 
cross-sections (<tnn) can be sorted out by 
noticing the sensitivity of Ebai on the system 
size, impact parameter as well as isospin de- 
gree of freedom of the reaction [6-10]. 
Experimentally, balance energy is observed 
for different systems ranging from 12 C + 12 C 



to 



197 



Au + 197 Au I6H 



15 



161 ] . The very ac- 
curate measurement of the E^ai in 197 Au+ 197 
Au 



16| has generated a renewed interest in 
the field. Unfortunately, these studies have 
not provided any significant contribution of 
isospin effects towards the balance energy. 
Later on, Pak et al. [8] demonstrated the 
isospin effect on the collective flow and bal- 
ance energy at central and peripheral ge- 
ometries. These findings were limited only 
for 58 Fe + 58 Fe and 58 Ni + 58 Ni systems. 
Theoretically, the disappearance of directed 
flow is studied using the Boltzmann Uehling 
Uehlenbeck (BUU) model. Q, H, 12, [l5| 



and Quantum Molecular Dynamics (QMD) 
model [3-5,10-21]. Different theoretical at- 
tempts considered either a stiff or soft equa- 
tion of state along with a variety of NN cross- 
sections. Very recently, Puri and co-workers 
flo| . conducted a very detailed analysis on 
the balance energy over entire periodic ta- 
ble with masses between 24 and 394. These 
study shed light on various aspect of nu- 
clear dynamics. Unfortunately, this study 
along with all other studies reported in lit- 
erature are limited to Central/Semi- Central 
Collisions only 10] . Following this work, the 
detailed analysis on the the semi-central and 
peripheral collision is performed by the Puri 
and co-workers in 2010 22|. All these studies 
indicated enhanced cross-section of 40-55mb 
with stiff equation of state to verify the bal- 
ance energy in intermediate energy heavy-ion 
collisions. All these studies were independent 
of isospin effects. The first study showing the 
isospin effects on the collective flow and bal- 
ance energy was reported by Li et al. [pj us- 
ing the isospin dependent Boltzmann Uehling 
Uehlenbeck (IBUU) model, where strong de- 
pendence of isospin effects was observed. In 



another contribution 



131 ] . they suggested the 



demand of reduced isospin dependent cross- 
section^ = 0.88o"7VAr) to better explain the 
experimental data. Chen et al. [ll| stud- 
ied the effect of isospin degree of freedom 
on the balance energy using isospin depen- 



2 



dent Quantum Molecular Dynamics (IQMD) 
model, which was an improved version of 
original QMD model Q, Q, 18] . The calcu- 



lated results were found to differ from the 
data at all colliding geometries. Recently, 
Gautam et al. [23j also studied the isospin 
effect on the balance energy by using IQMD 
model 24|, and also compared their find- 



ings with the other theoretical and exper- 
imental findings. They demanded to take 
care the Gaussian width(L) and cross-section 
(a = 0.88o"ATiv), while studying the isospin ef- 
fects in intermediate energy heavy-ion colli- 
sions. 

From the above, it is cleared that one is de- 
manding an enhanced constant cross-section 
[lo| in QMD model. On the other hand, it is 
also observed that the reduced isospin depen- 
dent cross-section is valid for the soft as well 



as for so 
state 13 



t momentum dependent equation of 



23j within IQMD or IBUU model. 



The systematic concept of enhanced and re- 
duced isospin dependent cross-sections in the 
presence of hard equation of state is missing 
throughout the literature. Moreover, the ef- 
fect of isospin dependent cross-sections are 
studied on the limited systems experimen- 
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23 



24). 



tally as well as theoretically [8|, [9 
We plan to study the effect of enhanced (30 
%of <Jnn) as well as reduced (30 % of &nn) 
isospin dependent cross-sections in the pres- 
ence of hard equation of state on the systems 



for which the experimental finding energy is 
predicted in the literature and then will com- 
pare the results with experimental findings. 
For this study, we will employ isospin depen- 
dent Quantum Molecular Dynamics(IQMD) 
model which is discussed in sec. II. The re- 
sults are discussed in sec. Ill, followed by 
conclusion in sec. IV. 



II. ISOSPIN-DEPENDENT QUAN- 
TUM MOLECULAR DYNAMICS 
(IQMD) MODEL 

The isospin-dependent quantum molecu- 
lar dynamics (IQMD) [24( model treats differ- 
ent charge states of nucleons, deltas and pi- 
ons explicitly 24| , as inherited from the VUU 
model [25)]. The IQMD model has been used 
successfully for the analysis of large number 
of observables from low to relativistic ener- 



gies 



231 ] . The isospin degree of freedom en- 



ters into the calculations via symmetry po- 



tentia 
tions 



cross-sections and Coulomb interac- 



251 ] . The details about the elastic and 



inelastic cross-sections for proton-proton and 
neutron-neutron collisions can be found in 
Ref. 24|. These cross-sections follow the 



data published by particle data group (PDG) 
for proton-neutron and proton-proton scat- 
tering j^J . In this model, baryons are repre- 
sented by Gaussian-shaped density distribu- 



3 



tions tions of motion: 

dfl d( H ) dpi d( H ) 

ffjtp t \ = 1 c -(?-fl(t)) 2 ^ C -(P-Pi(t)) 2 $_ dt dpi ' dt dri 



(2) 



(1) 

Where L is the Gaussian Width. As men- 
tioned in Ref. 



with 



= {T) + (V} 
24|, in IQMD the value of 2 



Gaussian width L depends on the size of 

i i j>i 

the system. This system size dependence xfj(f t ,p^,t)dfdr > dpdp f . (3) 

of L in IQMD has been introduced in or- 
der to obtain the maximum stability of the 
nucleonic density profile. Therefore, in the 



The baryon-baryon potential V % \ in the 
above relation, reads as: 



present study, by checking the stability, we T/^'iV _ r) — 4- 4- 4- 

\ J Skyrme 1 Yukawa 1 Coul 1 syn 

have taken the value from 0.5L to L. Its ear- / , . ... . n ( ' f* ' + f 

= ( ti^^ - r) + t 2 5(r^ - r)p 7 

lier version QMD has been very successful in 



explaining the multifragmentation [27j, tern- + £ 



exp(\f" — r\/fi) ZiZje 



2 



(Ir* — r|//i) |r* — r| 

prature and density [28|], flow |29|], multifrag- 1 / 

+t 6 — T;Ti5(fi -fj). (4) 



ments [2( and particle production 30(. 

Nucleons are initialized in a sphere with ra- Here Zi and Zj denote the charges of i th and 

dius R = 1.12A 1//3 fm, in accordance with the j th baryon, and Tj, T| are their respective 

liquid drop model. Each nucleon occupies a T3 components (i.e. 1/2 for protons and -1/2 

volume of h 3 , so that phase space is uniformly for neutrons). Meson potential consists of 

filled. The initial momenta are randomly Coulomb interactions only. The parameters 

chosen between and Fermi moment um(Pp). /i and t\, , t$ are adjusted to the real part 

The nucleons of target and projectile inter- of the nucleonic optical potential. For the 

act via two and three-body Skyrme forces, density dependence of nucleon optical poten- 

Yukawa potential, Coloumb interactions. In tial, standard Skyrme-type parametrization 

addition to the use of explicit charge states of is employed. The choice of equation of state 

all baryons and mesons a symmetry potential (or compressibility) is still controversial 

between protons and neutrons corresponding one. Many studies advocate softer matter, 

to the Bethe- Weizsacker mass formula has whereas, much more believe the matter to 



been included. be harder in nature 



25 



311 ] . For the present 



The hadrons propagate using Hamilton equa- analysis, a hard (H) equation of state, has 



been employed along with standard energy pendent cross-sections. We have attempted 
dependent cross-sections. Note that the to fit the reduced isospin dependent cross- 
relativistic effects are neglisible at these sections in the presence of stiff equation of 



enegies 



321 ] . state with experimental findings, as is per- 

formed in the literature with soft equation 
of state with and without momentum depen- 
dent interactions. 



There are two methods in the 



iterature used 



III. RESULTS AND DISCUSSION 

We study the directed flow using a stiff to find the balance energy [lflj. In the first 

equation of state along with enhanced and case ' tne balance energy is extracted from the 

reduced isospin dependent cross-sections (0= (-^M) P lots > where (P x /A) is to plotted as a 

0.7 to 1.3 <t nn ), by simulating various re- function of rapidity distribution Y c . m ./Y beam , 

actions. The time evolution of reaction is wn i° n * s gi ven as - 

follow upto 200 fm/c. This is the time at Y(i) = - In —^- (R) 

■ i r 2 E{i)- p S) 
which transverse m-plane flow saturates for 

n . t, ii r i • , n xi • where E(i) and P-Ai) are respectively, the to- 

hghter as well as tor heavier systems, rorthis w w 

study, the reactions of 12 C 6 + 12 C 6 (S = 0.4, tal energy and lon S itudinal momentum of i* 

L=0.5L) where L=8.66 fm\ ™Ne 10 + *Al n partide - Naturally the ener ^ at which this 

(6 = 04 L=0 5L) 40 Ar 18 + 45 Sc 21 (6 = fl ° W P asses tnrou S n zero is called balance en- 

4 L=0 5L) A0 Ar 18 + 51 V 23 (6 = 03 ergy ' The second method is to stud Y the inci - 

L=0 5L) m Kr m + 93 Nb 41 (6 = 04 dent ener Sy dependence of the directed trans- 

T n m \ 64 5s at- n ox verse in-plane flow (P* r ), which is defined as 

L=0.6L), Zn 30 + Ni 28 {b = 2jm, i — i r \ x n 

L=0.6L), 93 Nb 41 + 93 Nb 41 (6 = 0.3, ^ 

L=0.7L), 129 Xe 54 + ll8 Sn 50 (6 = 0- { P f) = \\~]sign{Y{i)}P x {i) (6) 

3/m, L=0.7L), 139 La 57 + 139 La 57 (6 = 0.3, A i 

L=0.8L), 197 Au 79 + 197 Au 79 (6 = 2.5 fm, where Y(i) is the rapidity distribution as dis- 

L=L) are simulated. The choice of impact cussed above and P x (i) is the transverse mo- 



parameter is guided by t 



re experimentally mentum of the i th particle in x-direction. 



extracted information j^sl, 3 The This (P^ vr ) is defined over entire rapidity re- 
above reactions were simulated between 45 gion and therefore expected to present an 
and 200 MeV/nucleon using the hard equa- easier way of measuring the in-plane flow 
tion of state along with different isospin de- rather than complicated (P x /A) plots. In the 
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FIG. 1: The averaged (P x /A) as function of the 
rapidity distribution. Here we display the result 
for Kr+Nb system at different incident energies 
and different isospin-dependent cross-sections. 

present study, we have tried to study the ef- 
fect of isospin dependent cross-sections on the 
flow or alternatively on the balance energy by 
using both of the parameter and then the de- 
tailed study is extended with later one. 

In Fig. 1, we display the change in 
the transverse momentum (P x /A) as a 
function of the rapidity distribution at 
different incident energies from 60 to 200 
MeV/nucleon for 86 Kr 3G + 93 iV6 41 systems. 



-2-101-2-1012 

Y /Y. 

cm. beam 



FIG. 2: As in [H but for different system and 
at particular cross-section a = 0.9cx/vAr- Diffrent 
panels are at diffrent incident energies. 

The different lines in the figure are showing 
the variation with different cross-section 
values. From the figure, we see that slope 
becomes less negative or more positive with 
increase in the incident energy. On the other 
hand, with reduction in the cross-section 
(o'nn), slope is getting more negative or less 
positive, while, becoming more positive with 
enhanced cross-section. This indicates that 
we see a change in the slope with incident 
energy and reduced isospin dependent cross- 
section. The figure is indicating two values 
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FIG. 3: Time evolution of (Pf) for different FIG. 4: Time evolution of (P£ ir ) for different 
systems at E=60 MeV/nucleon. The different systems at E=200MeV/nucleon. The different 
lines in figure representing the directed fow at lines in figure representing the directed fow at 



different cross-sections. 



different cross-sections. 



of balance energy i.e. E=80MeV/nucleon (at The slope is becoming more positive or less 

a = 0.7a nn) and around E=60MeV/nucleon negative with increase in the composite mass 

(at a = 0.9cttvjv) for 86 -K> 36 + 93 iV6 41 of system, indicating that lighter systems 

system. As the experimental balance energy remain in the environment of mean field 

for S6 Kr 3e + 93 iV6 41 is in the range of 55-60 compared to NN collisions at any given 

MeV/nucleon, So one is expecting to follow incident energy. The contribution of the 

the whole dynamics at a = 0.9<jnn- mean field the and collisions is discussed in 



Further, the detailed analysis of rapidity detail in Ref. 



10]. One also notice that a 



distribution of transverse momentum (P x /A) higher incident energy is needed in lighter 
for different systems at different energies cases to balance the attractive and repulsive 
with a = 0.9o"atjv is displayed in Fig. 2. forces. This energy is supposed to decrease 
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FIG. 5: Energy dependence of the directed nu- 
clear flow (Px ir ) for different systems. The lines 
have same meaning as that in fig 3 and fig 4. 

with increase in the system mass. Similar 
findings are also published by Puri and 
co-worker 10[ . Note that their study did not 



take isospin effects into consideration. 
Figures 3 and 4 are displaying the 
time evolution of the second parameter 
(P^ ir ) below (60MeV/nucleon) and above 
(200MeV/nucleon) the experimental balance 
energy, respectively. The results in the 
figure are displayed for five different systems 
and at reduced as well as enhanced isospin 
dependent cross-sections. The figures are 



FIG. 6: same as in fig 5, but for different sys- 
tems. 

indicating the similar scenerio with (-P* r ) as 
is depicted with (P x /A). Below the balance 
energy (in Fig. 3), the directed in-plane 
flow is negative during the initial phase of 
reaction for all the systems under considera- 
tion. This becomes positive at sufficient high 
incident energy say E=200 MeV/nucleon (in 
Fig. 4). These results shows that interaction 
among nucleons are attractive during the 
initial phase of the reaction, which turns out 
to be repulsive with increase in the incident 
energy. These interactions remain either 
attractive or repulsive throughout the time 
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FIG. 7: The shift in balance energy due to cross 
section as a function of combined mass of the 
system. 

evolution depends on the incident energy, 
isospin dependent cross-sections as well as 
composite mass of the system. It is clear 
from the figure that directed flow is becoming 
more positive or less negative with incident 
energy, isospin dependent cross-sections as 
well as with size of system. There is a sharp 
transition for each system from negative 
to positive directed flow at a particular 
cross-section. This particular transition 
is not possible for 12 C 6 + 12 C 6 system 
indicates the requirement of other variable 



FIG. 8: Balance energy as a function of com- 
bined mass of the system. The experimental 
points are represented with stars, QMD+40mb 
with crossed triangle and present with solid 
square. 

like momentum dependent interaction as 
well as enhancement in cross-section more 
then 30 %. If one compares the figs. 1-4, the 
same physics of balance energy elaborates 
with the (P x /A) as well as with (Pf r ). Out 
of these, as discussed earlier, (-P^ r ) is the 
simple and more useful quantity, because it 
is summed over entire rapidity distribution, 
that is why, (P^ ir ) is elaborated in detail for 
the further study. 
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FIG. 9: N/Z dependence of balance energy at 
a = 0.9cr jvjv The curve is parametrized with 
power law. 

To study the influence of reduced as well 
as enhanced cross-sections on directed 
flow (-P* r ) or alternatively on the balance 
energy, in Figs 5 and 6, incident energy 
dependence of directed flow is displayed 
for different systems. The different lines in 
figure represent the variation with different 
cross-sections. The studies with enhanced 
and reduced cross-sections are also avail- 
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The 



able in the literature 
experimental data are represented by stars. 
The directed flow goes from negative to 
positive value with increase in the incident 



energy. This is the general trend and is 
explained many times in the litrature by 
taking the concept of mean field and NN 
cross-sections. On the other hand, the 
role of different cross-sections is consistent 
through the present mass range. By finding 
the evidence of reduced cross-sections from 
Fig. 5, the results are displayed between 
(0.7-0.9 <x./v7v ) values in Fig. 6. The enhanced 
cross-section (1.3 ct^n) gives more positive 
value followed by the cugnon cross-section 
(ctnn) towards the reduced cross-section (0.7 
<tnn)- In other words, with increase in the 
cross-section value from (0.7 a nn-^-^nn), 
the directed flow is becoming more positive 
or less negative. This is due to the reason 
that with increase in the cross-section 
value, probability of reaction to take place 
increases that further results increase in 
the NN collisions and hence more positive 
value of the directed flow. This is resulting 
decrease in the balance energy. The balance 
energy is also found to decrease with increase 
in the composite mass of the system. This 
is due to dominance of Coulomb repulsion 
with an increase in the composite mass of 
system. Except for some lighter systems, the 
cross-section a = 0.9ctnn is found to explain 
the experimental balance energy nicely. 
Similar parametrization was also performed 
by Sood et al. [lO] within QMD model. 
They also found that enhanced cross-section 



10 



(cr = 40m6) can best explain the data, it goes on decreasing with system mass. It 

In contrary, calculation in IQMD model is also indicating theoretically that balance 

demand reduced value of cross-sections, energy is almost independent of the nucleon- 

The difference is due to the additional nucleon cross-section for the heavier system 

effect of isospin dependent cross-sections such as Au+Au, U+U etc. 

in IQMD model 241 ] . which were absent In Fig. 8, we display the energy of vanishing 

in QMD model. As in QMD model, the flow or balance energy (Ef, a i) as a function of 

strength of nn, pp, pp cross-section is taken composite mass of system that ranges from 

equal, while in IQMD, a np = 3a pp m 3a nn 40 v4r 18 + 45 Sc 2 i to 197 Au 7Q + 197 Au 79 . In this 

24j. Due to the different strength of np, figure, E hai is showed for the experimental 

pp, nn cross-section in IQMD, additional data (open stars), QMD+40mb (crossed tri- 

repulsion is produced compared to QMD angle) and IQMD+0.9 cjnn (solid squares), 

model. This addition repulsion will force the All the curves are fitted with power law 

directed flow to take earlier transition from of the form C(A TO t) t - The experimental 

negative to positive value and hence will data are fitted by r = —0.33 ± 0.06 The 

lower the balance energy in IQMD model balance energy is found to decrease with 

as compared to QMD for same cross-section the composite mass of the system, which is 

value. That is why, the balance energy that a well known trend discussed many times 



was obtained with QMD at a = AOmb are in literature [10J. The difference is in the 

at a = 0.9cr7VAr in IQMD model. This is r values obtained by different theoretical 

first ever parametrization of balance energy model. The BUU model report r between 

with hard equation of state in the presence —0.28 < r th < —0.32. In another study js] 

of reduced cross-sections with experimental again with BUU model r th = —0.41 ± 0.03. 

available balance energy. The present calculation depicts the r value 



By taking the Ref. 10j into account, which (—0.29 ± 0.06), which is close to the experi- 

depicts, that for heavier colliding nuclei E ba i mental r value (—0.33 ± 0.06) as compared 

is independent of the cross-section one is to QMD+40mb calculation having r value 

choosing, we have plotted in Fig. 7 the (—0.27 ± 0.17). In other words, the present 

AE ba i = {E ba i) , 7aNN - {E ba i) , 9(TNN with IQMD model with a stiff equation of state 

composite mass of system. Our findings are along with a = 0.9cr tvv can explain the 



also supporting the findings of Ref. 



12j. data much better than any other theoretical 



AEh a i is maximum for lighter systems and calculations. The a = 0.9<jnn explains the 
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data for all nuclei, except for some lighter A large number of reactions were studied 
nuclei. The lighter nuclei, when checked having mass range from 24 to 394, where 
out, demand for an enhanced cross-sections experimental balance energy is available. 



3, 



33| along with momentum dependent Our calculation with stiff equation of state 



interactions 



27| . Our calculations about and reduced cross-section (a = 0.9<jnn) are 



the strength of reduced NN cross-section is in good agreement with the experimental 
in agreement with earlier calculation, where findings, except for 12 Cq + 12 Cq. The depen- 
disappearance of transverse in-plane flow dence of isospin dependent cross-sections get 



23[ as well as elliptical flow is parametrized weakens with increase in the size of system, 

with experimental data [13]. The balance energy is parametrized with 

We have also tried to fit the balance energy N/Z ratio in terms of power law, which is 

in terms of other parameter such as the to be quite similar with the parametrization 

neutron to proton ratio of colliding nuclei, of composite mass of system, but the r 

This attempt is shown in Fig. 9, where values are different in both of the cases. One 

balance energy is plotted as a function of could try the balance energy prediction with 

N/Z. The Eb a i is parametrized with power enhanced isospin dependent cross-section 

law of the form (N/Z) T . The r value in in the presence of momentum dependent 

N/Z dependence is —2.39 ± 0.40, while in interaction for 12 Cq + 12 which is earlier 



10| by using QMD 



Atot dependence in Fig. 8 is —0.29 ± 0.06. studied by Sood et al. 
The t value in this case is larger compared in the presence of momentum dependent 
to the mass dependence. The difference in interactions, 
the slopes may be due to different charge to 
mass ratio in heavier colliding nuclei. 



IV. CONCLUSION 

By using the IQMD model, we have 
studied the effect of reduced as well as 
enhanced isospin dependent cross-sections 
on the directed flow and balance energy. 
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